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Abstract 

The urgency of predictions in large-i region at LHC stimulated us to present a unified model of small and high t 
scattering at high energies. Our model is based upon a safe theoretical ground: analyticity, unitarity, Regge behavior, 
gluon exchange and saturation of bounds established in axiomatic quantum field theory. We make precise predictions 
for the behavior of the differential cross sections at high t, the evolution of the dip-shoulder structure localized in the 
region 0.5 < \t\ < 0.8 GeV 2 and the radical violation of the exponential behavior of the first diffraction cone at small t. 
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R.F. Avila, P. Gauron and B. Nicolescu (AGN) 
recently published a good fit of forward and non- 
forward data in pp and pp scattering, leading to 
intriguing predictions at RHIC and LHC p] . 

The AGN model was defined to be valid only 
at moderate-i values (0 < \t\ < 2.6 GeV 2 ). We 
stress that this model is based upon a safe theoret- 
ical ground: analyticity, unitarity, Regge behavior 
and saturation of bounds established in axiomatic 
quantum field theory. 

In the present paper we generalize this model by 
extending its validity till high-i values (0 < \t\ < 
16 GeV 2 ) and we present our predictions at RHIC 
and LHC. Such predictions are urgently needed in 
view of the near start of experiments at LHC. 

As data we used the dataset recently proposed 
by J.R. Cudell, A. Lengyel and E. Martynov 
(CLM), who built a coherent dataset of all exist- 
ing data for 5 < -/s < 1800 GeV and < \t\ < 16 



GeV 2 [2] . The raw data are, of course, contained in 
the Durham Database [Jj, but CLM introduced a 
detailed study of the systematic errors and gath- 
ered into a common format more than 260 subsets 
of data from more than 80 experimental papers. 
The CLM can be considered as a standard dataset, 
which is available online [I] . 

Let us first define the even-under-crossing and 
the odd- under-crossing amplitudes F±(s,t): 



F± = -(F pp (s,t)±F Pp (s,t)) , (1) 



which are normalized so that 
1 



0t(s) 



y / s(s-4m 2 ) 
ReF(s,t = a) 
ImF(s,t = 0)' 



ImF(s,t = 0), 



(2) 
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^ ( t) ^ 
dt ' 16ngs(s — 4m 2 ) 

g = 0.3893797 mb- GeV" 2 , 
m p being the mass of the proton. 



\F(s,t)f 
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In the AGN model" 2 ! F + (s,t) is defined as a 
superposition of the following contributions: the 
Froissaron Fp(s,t), the Pomeron pole F p (s,t), 
the Pomeron-Pomeron Regge cut F pp (s,t), the 
secondary / 2 — a 2 reggeon F+(s, t) and the 
reggeon-Pomeron cut F pp (s,t): 

F+{s, t) = F?(s, t) + F p (s, t) + F pp (s, t)+ 
F p (s,t)+F pp (s,t) , 



(4) 



where 



Ff (s,t) = is { Hx In 2 § 2Jl ^ K ^ exp(b +1 t) 



K+T 



-if 2 In s Jo {K+t ) In s exp(& + 2^) 
-H 3 [J (K+f) - K+tJ^K+t)] exp(6 +3 i)} 



F p (s,t) = C P exp(bpt)* + (a P {t)) 

(s/s 



,a P (t) 



(5) 



(6) 



F pp (s,t) = ^ILeMbppt)<i> + (app(t)) 
In s 

( S / So ) QppW , (7) 
F p (s,t) = C+ exp(b+t)<f> + (a+(t)) 

( S /s r>\ (8) 

F* p (s,t) = Z^E exp ( & + i3 t)$ + (a+ p (t)) 

(s/soT^W, (9) 



Ins 



1 



.5 = (s/s ) exp ( --»7T ) , s = lGeV^, (10) 



f = y/-t/to1n.S, t = lGeV 2 , (11) 
*+(«(<)) = *sin ~ cos (|a(*)) , (12) 



ap(t) = 1 + apt , 
Oipp(t) = 1 + a'pt/2 , 
a+(i)=a+(0)+a+'t, 



(13) 
(14) 

(15) 



4pW = «h(o) + -rrV* > ( 16 ) 



a p + a R 



Hi, H 2 , H 3 , K + , b+i, b +2 , 6+3, Cp, bp, Cpp, 
bpp, Cp, frpp, Cpp, b\ p , a' P , aj(0), ap; being 
constants. Jq and Ji are Bessel functions. 



2 Notation for the signature factors differs from the one 
of the original AGN model. 



In its turn, the F-(s,t) is defined, in the AGN 
model, as a superposition of the Maximal Odd- 
eron Contribution F^ Io {s,t), the Odderon pole 
F°_(s,t), the Odderon-Pomeron cut F® p (s,t), 
the secondary p — lo reggeon F p (s,t) and the 
reggeon-Pomeron cut F pp (s, t): 



F_(s,t) = F_ + F_(s, t) + F_ (s, t)+ 



F p {s,t)+F RP (s,t) 



(17) 



where 



F?° (s,t)=s{ 0l In 2 s^f) 



exp(6_i£) 



+O2 In s cos(/XT_f ) exp(6_ 2 t) 
+O3 exp(fe_ 3 t) } . 

F°(s,t) = C7 exp(M)*-(ao(*)) 

(»/*o) 

F OP (s,i) = ^exp(6 pi)$_(a p(t)) 



oo(t) 



Ins 



( s / s o) 



aop(t) 



(18) 



(19) 



(20) 



F*{s,t) = -Cpexp(&pi)$_(ap(i)) 



In s 



(s/s ) Q «W, (21) 

exp(&ppt)$_(a^ P (t)) 

(s/s ) Q «- (t) , (22) 
$_(a(f)) = icos(^a(*)) + sin (23) 

ao(t) = l + a' t ) (24) 



aop(i) = 1 + — 



'p"o 



a P + a. 



(25) 
(26) 



«5p(*)=«5(0) + -rT^T*. (27) 
a P + Q fl 



Oi, 2 , 3 , K_, 6 1 , b 2 , b 3 , C , b Q , C p, 
bop, Cp, b~, Cpp, bpp, a' , a R (0), ap being 
constants. 

The AGN amplitudes defined through eqs. (JTJ- 
(|27| , lead to an excellent description of the stan- 
dard dataset: 



X 2 / do f =1.16 



(28) 



for a total of 1966 experimental points. The ap and 
p data (238 experimental points) correspond to an 
energy range 5 GeV < yfs < 1.8 TeV. The da/dt 
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Table 1 



Table 2 



The values of parameters in F+(s,t) 



Name Dimension Value Error 



Hi 
n 1 




GeV 2 


0.22 


0.002 




mb 


GeV 2 


-0.016 0.001 


Hi 


mb 


GeV 2 


27.43 


0.30 


K 






0.20 


0.001 


c, 


mb 


GeV 2 


2.65 


0.16 


Cpp 


mb 


GeV 2 


-40.04 


2.18 


C+ 


mb 


GeV 2 


77.61 


1.71 


C+ 


mb 


GeV 2 


-23.73 


0.70 


OL py J ) 






0.67 


0.01 


a+' 


GeV" 2 


0.80 


0.01 


a p 


GeV" 2 


0.41 


0.01 




GeV - 2 


3.06 


0.03 


6+2 


GeV" 2 


0.66 


0.01 


6+3 


GeV" 2 


4.57 


0.03 


bp 


GeV" 2 


0.31 


0.07 




GeV" 2 


9.18 


0.44 


b+ n 


GeV" 2 


4.80 


0.14 


b RP 


GeV" 2 


0.46 


0.02 


N+ 


mb 


GeV 2 


-9.42 


0.33 


t+ 


GeV 2 


0.45 


0.01 



data (1728 experimental points) correspond to an 
energy range 9.3 GeV < yfs < 1.8 TeV and to a 
i-range 0.1 < \t\ < 2.6 GeV 2 . The presence of the 
Maximal Odderon is, of course, crucial in order to 
describe the experimental difference between the 
pp and pp differential cross-sections at y/s — 5.8 
GeV in the dip-shoulder region 1.1 < \t\ < 1.5 
GeV 2 0. 

Now, we generalize the AGN model by introduc- 
ing the following three ingredients: 
1. Contributions N±(s, t) to F±(s,t) behaving like 
t~ A at high t: 

N+(s, t) = isN + In — L— , (29) 



N-(s,t) = siV_ Ins 



(l-i/*+) 6 ' 
-t 



(30) 



(I-*/*-) 5 ' 

These additional terms add just 4 free parameters 
more as compared with the original AGN model. 
Similar forms as those of Eqs. (|2"9"|) - (f30"|) were used 
inRef. [6]. 

The experimental t~ A behavior at large t is well 
known. 

The theoretical motivation for N-(s,t) was 
given longtime ago by Donnachie and Landshoff: 



Name Dimension Value 


Error 


Oi 


mb 


GeV 2 


0.002 00.001 


o 2 


mb 


GeV 2 


0.12 


0.01 


o 3 


mb 


GeV 2 


-2.13 


0.08 


K_ 






0.011 


0.002 


Co 


mb 


GeV 2 


0.88 


0.02 


Cop 


mb 


GeV 2 


0.061 


0.005 


R 


mb 


GeV 2 


38.00 


0.37 


RP 


mb 


GeV 2 


-4753 


101 


q"(0) 
R \ 1 






0.53 


0.01 


R 


GeV" 2 


0.80 


0.01 


(J 


GeV" 2 


0.10 


0.01 


b-i 


GeV" 2 


1.07 


0.02 


6-2 


GeV" 2 


2.87 


0.06 


6-3 


GeV" 2 


3.09 


0.04 


bo 


GeV" 2 


2.29 


0.02 


boP 


GeV" 2 


0.18 


0.01 


IX 


GeV" 2 


1.42 


0.05 


b RP 


GeV" 2 


8.44 


0.77 


N- 


mb 


GeV 2 


2.16 


0.20 


t- 


GeV 2 


0.61 


0.01 


A-MO 


GeV" 2 


7.32 


0.32 


Ao 


GeV" 2 


8.87 


0.41 


t c 


GeV 2 


-0.16 


0.004 



the Odderon 3-gluon exchange, with elementary 
gluons [7]. The motivation for N+(s, t) is more ob- 
scure, but it could be correlated with the C = + 
part of the exchange of 3 elementary gluons. 
2. A crossover factor Z-(t) in the F_(s,t) ampli- 
tude 

tanh(l - t/t c ) 



z-(t) = 



(31) 



tanh(l) 

This factor introduce just one more free parameter 
as compared with the original AGN model. 

The phenomenological motivation for such a 
Z-(t) factor is obvious: the pp and pp cross at 
very small t (\t\ ~ 0.16 GeV 2 ) and the respective 
crossover point is practically independent on s. 
3. Two linear functions (1 + Awot) an d (1 + Aot) 
which multiply the Maximal Odderon and the 
Odderon pole contributions, respectively. These 
factors have also a phenomenological justification: 
they allow us to better describe the smallness of 
the Odderon forward couplings at present ener- 
gies. These linear functions introduce just 2 more 
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Fig. 1. Total cross-sections. Our predictions at RHIC and 
LHC are also shown. 




10 V(GeV) 1 ° 3 



Fig. 2. pp elastic cross-sections. Our predictions at RHIC 
and LHC are also shown. 

free parameters as compared with the original 
AGN model. 

In fact we considered a big number of fits by tak- 
ing into account just one additional ingredient or 
a combination of two additional ingredients. The 
net result is that, for an excellent fit of the stan- 
dard dataset we need all the 3 ingredients. This 
result is not trivial, because we add a ridiculously 
small number of free parameters (7 free parame- 
ters) as compared with the total of 36 free param- 
eters of the AGN model. Moreover, we consider- 
ably increase the number of fitted data points: from 
1966 to 2892. Therefore, we are convinced that the 
3 ingredients represent a safe theoretical general- 
ization of the AGN model at high t. 

The best value of the free parameters involved in 
the amplitudes defined in Eqs. l|4"])-([2"7]l and 



n"j) with the additional 3 ingredients are obtained 
through a x 2 MINUIT minimization and they are 
shown in Tables 1 and 2. 

The x 2 /dof is surprisingly small: 

X 2 /dof = 1.23, (32) 

by taking into account the huge number of exper- 
imental data which are fitted (2892 points) in a 
huge range of energies 

5 < Vs < 1800 TeV (33) 

and in a huge range of rl 

t = and 0.1 < \t\ < 16 GeV 2 (34) 

The fit presented in this paper is certainly the best 
existing model for describing the data. We present 
here all our analytic forms and all the values of 
the parameters so any author can compare his /her 
own model with our model. 

The quality of the fit is shown in Fig. 1 (total 
cross-sections a tot (s)), Fig. 2 (elastic cross-sections 
er e /(s)), Fig. 3 (p-parameter p{s)), Fig. 4 (pp dif- 
ferential cross-sections), and Fig. 5 (pp differential 
cross-sections). In Figs. 1-3 we also show our pre- 
dictions for <7 tot (s), <J e i(s) and p(s) at RHIC and 
LHC energies. 

Our predictions for the pp differential cross- 
sections at RHIC and LHC energies are shown in 
Fig. 6. We predict (see Fig. 6) a dip - shoulder 
structure persisting at RHIC and LHC energies, 
namely: a dip at \t\ ~ 0.7 GeV 2 for y/s = 500 
GeV, which moves to \t\ ~ 0.5 GeV 2 for x/s = 14 
TeV, while the shoulder present at \t\ ~ 1 GeV 2 
for y/s = 500 GeV moves to \t\ ~ 0.8 GeV 2 at 
y/s = 14 TeV. 

In Fig. 7, we plot the slope B pp (s, t) of the differ- 
ential cross-sections at RHIC and LHC energies, 
as a function of s and t. A very interesting effect 
can be contemplated in Fig. 7: at RHIC energies 
B(s, t) is constant at fixed s in a quite large range 
of t but at LHC energies the constancy of B(s, t) 
is totally lost: the exponential behavior oida/dt is 
no more valid. 

All these predictions can be verified at RHIC 
and LHC energies. In particular, the TOTEM ex- 
periment will explore the large \t\ pp elastic scat- 
tering till 8 GeV 2 [5j. Also experiments at RHIC 
will be very helpful. The verification of our predic- 
tions can help to clarify the dynamics involved in 
small and high t scattering. 



3 The range in Eq. {34} is justified in Ref. We exclude 
244 points, most of them at low energies, from the dataset 
in Ref. [2] (corresponding to 5 pp sets and 2 pp sets), 
because they are in conflict with all the other data. 
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pp 



Fig. 3. p-parameter. Our predictions for pp scattering at 
RHIC and LHC are also shown. 



la (GeV) = 6.15 
9.78 
13.76 
16.82 
19.42 

23.50 (27.42) 
53.00 




Fig. 4. pp differential cross-sections in the energy range 
6.15 < y/s < 53 GeV. Data and theoretical values are 
multiplied by 10 — 2 (" — !), where n is the number of curve 
(and corresponding dataset) starting from the top. 



s (GeV) = 7.00 
9.78 
13.76 
53.0 
546.0 
630.0 
1800.0 




I .... I .'TTh 
12 13 14 15 16 



Fig. 5. pp differential cross-sections in the energy range 
7 GeV < y/s < 1.8 TeV. Data and theoretical values are 
multiplied by 10 — 2 ( n_1 ) ) where n is the number of curve 
(and corresponding dataset) starting from the top. 




Fig. 6. Evolution of pp differential cross-sections from RH1 C 
to LHC energies. 
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